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Abstract 



We examine variations in the equation of state of the cosmic string por- 
tion of the cosmological fluid which lead to perturbations of the background 
matter density. These fluctuations in the equation of state are due to varia- 
tions in the local density of cosmic string loops and gravitational radiation. 
Constructing a crude model of the distribution of entropy perturbations, we 
obtain the resulting fluctuation spectrum using a gauge-invariant formalism. 
We compute the resulting cosmic microwave background anisotropy, and esti- 
mate the effect of these perturbations on the cosmic string structure formation 
scenario. 
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I. INTRODUCTION 



Cosmic strings are topological defects which may have formed in a Grand Unified Theory 
era phase transition. Cosmic strings may play an important role in cosmology inducing 
temperature anisotropics in the cosmic microwave background (CMB), and seeding pertur- 
bations for the formation of large scale structure (LSS). 

In past work, Veeraraghavan and Stebbins [0] developed a formalism for the analysis of 
linearized perturbations to the cosmological fluid induced by cosmic strings. These tech- 
niques have been applied to the simulation of CMB temperature anisotropy 0, and to 
the analysis of string perturbations in cold dark matter for the formation of LSS In 
both applications, the effects of the strings on the linear density field were computed by 
convolving the string stress-energy tensor with the appropriate Green's function. As well, 
density perturbations, anti-correlated with the cosmic strings, were included to compensate 
for perturbations induced by strings at formation. 

We propose to compute the effects of an additional source of cosmic string-generated per- 
turbations, not explicitly included in previous analyses. These are entropy perturbations, 
caused by the inhomogeneous shifting of energy between cosmic string loops and gravita- 
tional radiation emitted by loops. The spatial variation in the equation of state of the total 
cosmological fluid is due to statistical fluctuations in the density of long cosmic strings, from 
which the distribution of small loops and gravitational radiation follow. Below a coherence 
scale ^, the equation of state of the loops and radiation is uniform. The loops and radiation 
act together as a dilute gas, shifting the expansion rate within a volume Above ^, the 
equation of state varies as rms statistical fluctuations, oscillating with wavelengths longer 
than the coherence scale. 

The reason why these effects have not been explicitly treated in past work are the fol- 
lowing. First, no simulation has included the gravitational radiation emitted by loops and 
strings in the cosmic string stress-energy tensor, for the purposes of CMB or LSS calcu- 
lations. The effect of gravitational wave, tensor perturbations on the scalar, gravitational 
potential is a second order effect. However, variations in the equation of state, as considered 
in this paper, result in a first order effect in linearized gravity. Second, the perturbations 
resulting from the variations in the equation of state, although not genuine fluctuations in 
the total density, have not been modeled by the compensations employed to date 0. The 
compensations included in CMB and LSS analyses have simply accounted for the initial 
distribution of strings and loops, not the subsequent shift in the form of matter, from loops 
and long strings to gravitational radiation. 

Our goal in this paper is to estimate the importance of entropy perturbations for the 
formation of large scale structures and CMB anisotropy, prior to carrying out a more detailed 
investigation (some preliminary work has been carried out in 0). In the following, we present 
a simple analysis of the effects of these perturbations from cosmic strings. In section |l| we 
present a model of the variations in the cosmic string equation of state. In section |IT1| we give 
explicit solutions for the gauge-invariant potential, characterizing the density fluctuations 
induced by the entropy perturbations. We also compute the power spectrum, and comment 
on the effects on large scale structure formation. In section |I^ we estimate the effect of 
these perturbations on the CMB. We conclude in section [V]. 
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II. MODEL OF ENTROPY PERTURBATIONS 



Entropy perturbations to the cosmological fluid arise when there is an inhomogeneous 
shift in energy from one form of matter to another. Consider a perturbing component of 
the fluid 5p with pressure 5p. The resulting perturbation to the entropy density is 

5S = T-\5p-cl5p) (2.1) 

where Cg is the speed of sound in the background medium at temperature T. Hence, the 
presence of an additional form of matter for which 5p ^ c^Sp, that is the perturbing matter 
does not obey the equation of state of the background fluid, leads to a perturbation of the 
entropy density. In the cosmic string scenario, the long cosmic strings, loops, and emitted 
gravitational radiation, 

Sp = Poo + Ploops + Pgr (2.2) 

serve as the perturbing component of the cosmological fluid, driving the entropy perturba- 
tions in ( |2.1|) . More formally, we write 

T6S = -Q^.g^" - (1 + c^)e^,M'^M^ (2.3) 

where Qfj^ij{ri,x) is the cosmic string stress-energy tensor, including the gravitational radia- 
tion, and is a unit time-like 4- vector orthogonal to the time-constant spatial hypersurfaces. 



Given the source (|2.3| ), we may carry out a procedure similar to that in convolve the 

source with the appropriate Green's function to obtain the gravitational potential. While it 
may be straightforward to adapt a cosmic string simulation to this task, in the present work 
we will construct an analytic model describing the cosmic string entropy perturbations. 

Treating the cosmic string network as a dilute gas, we may compute the pressure con- 
tributed by the strings. In the dust-dominated era, only the gravitational radiation con- 
tributes to the entropy perturbations, 5S in equation (|2.1|). Loops evolve as collisionless 
dust, so do not contribute. Hence, it is the shifting of energy from cosmic strings into 
gravitational radiation that drives the entropy perturbations. 

Based on a crude model of network evolution in which the strings rapidly approach 
a scaling solution, the energy injected into gravitational radiation, in the dust-dominated 
era, in a time interval given by the coherence time scale ^ is a constant fraction, JF, of 
the background energy density. Because the gravitational radiation is free-streamed at the 
speed of light, we estimate that the coherence scale is given by the Hubble length, ^ = 7i~^ . 
Hence, decomposing the perturbing energy density field in terms of an amplitude and spatial 
distribution, we write 



Here, t] is the conformal time, po is the background fluid energy density, JF is the rms ampli- 
tude of the statistical fluctuations in the energy density as a fraction of the background, and 
/(?7, x) is a normalized distribution describing the spatial variation of the entropy perturba- 
tions. In the dust-dominated era, J-'gr ~ AT^AGpil — 2(t>^)) ~ 0{\)Gp. In this expression 
A, defined by 
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A^^{{p^-{p^))l)e^^-\ 



(2.5) 



is the amphtude of the statistical fluctuations in the long string density poo, averaged on 
a length scale ^, where t is physical time. The squared, average velocity along the string 
is (f^) = 0.37 in the dust-dominated era 0. We caution that the decomposition ( |2.4| ) is 
a simplification. A more realistic analysis would use the stress-energy tensor of the cosmic 
strings, as determined by a numerical simulation, to compute 5S. 

The function /(r^, x) describes how the perturbations to the equation of state vary with 
position and time. We decompose this function in terms of fourier harmonics 

/(r/,x) = (27r)-3|d3^e^fc^-/(^,fc) 

{f{r^,k)f{i,k')) = {2iTf6{k-k')\f{r^,k)f{i,k)\, (2.6) 

where /(?], k) = |/(?7, fc)|e*^fc. The fluctuations in the density of long strings lead to statisti- 
cally homogeneous, isotropic variations in the density of loops and gravitational radiation. 
Because the long cosmic strings make a statistically random walk on scales larger than the 
string correlation length, the variations in the equation of state, and hence the phases 6j^, 
are randomly distributed. The product of the fourier coefficients, \f{ri,k)f{rj',k)\, plays a 
role similar to that of the structure function used by Albrecht & Stebbins to describe 
the matter perturbations directly induced by cosmic strings and loops. We normalize this 
distribution by smoothing the variance over a time scale ^: 

aj ^ (/(^, x)f{v\ x)) = 1. (2.7) 

The reason for smoothing in time is that the cosmological fluid can only react to changes 
in the local equation of state due to the cosmic string loops and gravitational radiation on 
time scales greater than the coherence time scale. Hence we need only restrict the variance 
of the distribution / on the time scales which will lead to entropy perturbations. 

We expect that the distribution does not vary on length scales smaller than the time ^, so 
that f{ri, k > ^^^) = 0. As a general distribution we write / as a power law in wavenumber 
k, on scales above ^: 

\f{v,k)\ = 27rv'(a + 5/2)(a + 3/2)e^/2(A;0"e(l - k^ (2.8) 

where ^/rj = constant. The numerical coefficients and dependence on ^ follow as a result 
of the normalization, in equation ( p.7| ). The step function is defined as 0(x) = 1 for a; > 
and Q{x) = for x = 0. The parameter a determines how scales above the coherence scale 
^ contribute to the distribution. As a — > oo 

jim |/(r/,fc)| = 27re=^/25(l-A;0 (2.9) 

so that the distribution is characterized by a single, comoving scale. For finite, decreasing 
values of a, the range of comoving scales contributing to / increases. In the limit a — > 0, 
the distribution has support on all wavelengths above the comoving scale ^. Because the 
network of long strings is distributed as a random walk on large scales, we expect the entropy 
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perturbations to similarly have support on scales above ^. Hence, we expect a to be 
more physically reasonable than a — >■ oo. However, in the following calculations we will take 
a — i> oo in order to simplify the calculations, noting that for other values of a, the results 
are not strongly affected. Having specified the properties of the entropy perturbations, we 
may proceed to determine the resulting fluctuation spectrum. 



III. EVOLUTION OF PERTURBATIONS 

To describe the perturbations of the background cosmological fluid, we use a gauge- 
invariant formalism [0,^. We must supply an expression for the density perturbation relative 
to the hypersurface which represents the local rest frame of the cosmological fluid; this pre- 
scription gives as close as possible a "Newtonian" time slicing. For an entropy perturbation, 
that is a perturbation to the equation of state, we find the evolution of the gauge-invariant 
potential $ is given by 

$" + 3(1 + cl)n^' + ^IJi! + (1 + "icDU^]^ - c^V2$ = ^TxGa^[bv - c^Sp). (3.1) 

Here we work with conformal time rj, where ' = d/dt], and Ti = a' /a. The coordinates [r], x) 
are dimensionless, as the expansion scale factor a carries units of length. 

We will also compute the power spectrum of the gauge-invariant potential generated by 
the entropy perturbations in the dust-dominated era. We define the correlation function of 
the potential to be the second moment of $: 

^V, r) = mv, xWv, x')) =471 f k^dk ^^^^^$(7/, k) (3.2) 

J kr 

Here r = \x — x'\ is a conformal distance. Recall that the density contrast 5 due to the 
gauge-invariant potential is given by 

5 = h-L-^ (V2$ - 3H$' - m'^^) . (3.3) 

For modes well within the horizon, k ^ Ti, the density contrast is dominated by the gradient 
term. Thus, we have the relationship 

Ak'^ 

n = ^Pc, for k:^n (3.4) 

relating the power spectrum in the potential to the power in the density contrast. 

In the dust-dominated era, the scale factor behaves as a oc rj"^ and the speed of sound 
for the cold, collisionless dust vanishes so that the evolution equation for the gauge-invariant 
potential simplifies greatly. We may solve the differential equation as 



$(77, f) = AttG r dr]'r]'-^ T dr]" r]"^a\r]"){^5p{r]" ,x) - c^SpiT]" ,x)) 



V rv 
drj'ri'-^ / 

+ {r] - r],g)<^'{r],^, x) + <i>{r],^, x). (3.5) 

Hence, given the boundary terms $ and $' at T]eq, at the onset of the dust-dominated era, 
it is straight forward to obtain the potential $ at any later time. 
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The source of the entropy perturbations is the gravitational radiation emitted by os- 
cillating cosmic string loops. Using (|2.4| ), we may obtain the dust-era solution for $. We 
will make the simplification that the boundary terms are unimportant. This is reasonable 
provided we are interested in length scales due to perturbations which enter the horizon 
after radiation- matter equality. Hence, we find 

AJ-"^ rv rv' ~ 2 

^*(^,A;) = ^ / d^'Wr'l dv"{vrfiv'\k) . (3.6) 

In the case a ^ oc, the power is 



2T'^ 

Jim P$(r/, k) = -^k-^[^ - x-^fe{x - 1)6(1 - x^y). (3.7) 

Here, x = fc^. For other values of a the power spectrum also behaves as "P^ oc k~^ for modes 
well within the horizon, generating a scale-invariant isocurvature spectrum. Fluctuations 
generated during the dust-dominated era will have wavelengths from ~ 10/i~^Mpc, the 
horizon scale at radiation-matter equality, up to the present horizon scale. These fluctuations 
will contribute to the large angle CMB anisotropy and the formation of the largest structures. 

We may compare the power in entropy-generated density fluctuations with the power 
due to the perturbations induced directly by the cosmic strings. Referring to ( p.4| ), we find 

ATrk'Vsik) = '^J'^rikhy (3.8) 

where k has units /i^/Mpc and h is the Hubble parameter normalized to Hq = lOOkm/s/Mpc. 
Compare this to the power spectrum due to cosmic strings in CDM (equation 8 of Q): 

Ank^Vsik) = 47c{6.8)\kh)\Gfif. (3.9) 

Here we have used the I-model of Albrecht & Stebbins |^ motivated by the string simu- 
lation results of Bennett & Bouchet and Allen & Shellard 0, on length scales well above 
the 8h~^Mpc normalization scale. Constructing a ratio of the LSS power due to the pertur- 
bations induced directly in the cosmological fluid by cosmic strings, to those generated by 
entropy perturbations, we obtain 

Vdirect/Ventropy ~ 50{Gfl)^J^-^ = 3.2A-\ (3.10) 

We see that for J^gr ~ 7Gfi or, using the estimate for J^gr preceding equation (P3|), for 
A ~ 1.8 the power for large scale structures are comparable. As a — 0, such that the 
entropy perturbations have support on a wide range of length scales, the power due to 
entropy perturbations increases slightly, and the powers are comparable for A ~ 1.4. The 
reason for the increase is that entropy perturbations generated at a range of times, rather 
than at a single time as with the a — >■ oo model, now contribute to the power in a single 
mode, on these scales. Numerical simulations indicate that the mean energy density in long 
strings during the dust-dominated era is poo = (4 ± l)/i/t^ 0, obtained by averaging over 
multiple realizations of the string simulation. If the quoted uncertainty in poo measures the 
amplitude of the statistical fluctuations, averaged on a scale ^, then A ~ 1. Modeling the 
long string network as a random walk on large scales, we expect A ~ 0.6 0. Hence, we 
expect A ~ 0.6 — 1 to be a reasonable range for the amplitude of entropy perturbations. We 
caution that these results are intended as a preliminary estimate; a more reliable calculation 
will use the string simulations directly to compute the effects of entropy perturbations. 
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IV. CMB ANISOTROPY 



The CMB anisotropies induced by the entropy-generated perturbations may be computed 
from the gauge-invariant potential. The dominant contribution on large angular scales, as 
is common for isocurvature perturbations, is due to the time variation of the potential $, 
integrated along the line of sight to the surface of last scattering: 

AT po 

— {n)=2 dri<!?'{ri,\rio-ri\n). (4.1) 

Here, n is a unit-vector on the celestial sphere. The correlation function is expressed in the 
usual form, 

AT AT °^ 2/ + 1 

m ^ {—{n) — {n')) = j:^C,P,{cose). (4.2) 

The amplitude of the multipole moments is given by C^, and h ■ n' = cosO. Evaluating 
the large-angle CMB anisotropy due to dust-era entropy perturbations only, in the a — cxd 
case, we find i^d ~ GJ^^^ for 10 < i < 40. Recall that J^gr ~ AAGfi, so that for A ~ 2.0 
these multipole moments are comparable to the f^C^ ~ 400 (G/i)^ perturbations measured 
in the numerical simulations [0]. This is roughly consistent with the result for y4(~ 1.8) 
found in the previous section, in equation ( p.lQ ), in order that the LSS power due to the 
entropy perturbations is comparable to the power due to the perturbations induced directly 
by cosmic strings. 

The effects of the entropy-generated perturbations on small-scale CMB anisotropy are 
estimated using the analytic techniques developed by Hu & Sugiyama 0. 

6i{r]o, k) = (Ooiriis, k) - ^{r]is, k)^ji{k\r]o - r]is\) 

+9,ir]is,k){2i+l)-'{ije-iik\7]o - Vis\) - {i + l)3i+i{k\Vo - Visl)) 

-2 rd7]^'{r],k)Uk\r]o-v\) (4.3) 

In the above equations (see equations 12 and 14 of 0), Oq and 6i are the monopole and dipole 
moments of the perturbations to the baryon-photon fluid. The final term, the integrated 
Sachs- Wolfe effect, dominates on large angular scales, as with isocurvature models. Modeling 
only the dust-era entropy perturbations, the resulting contribution to the CMB anisotropy 
due to monopole and dipole oscillations of the baryon-photon fluid is small. However, if 
we artificially lengthen the duration of the pre-recombination, matter-dominated epoch as a 
way of extending the validity of our model to smaller scales and earlier times, we find that 
significant anisotropy may result on small angular scales. For the a ^ oo model with a 
coherence scale ^ = H~^, a series of peaks develops, beginning near i > 200, before diffusion 
damps out the perturbations. Smaller values of the coherence scale lead to a decrease in the 
the peak amplitude and a shift in the peak location to higher i values. It is important to 
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note that the location of the peak in the spectrum is sensitive to the value of the comoving 
coherence scale ^ for the entropy perturbations. Magueijo et al |10] have carried out a 
similar calculation of the small-angle anisotropics using the Albrecht-Stebbins LSS model as 
an effective source for fluctuations in the baryon-photon fluid. There, they found that the 
location of the peak shifts to larger values of £ as the coherence length scale decreases. This 
behavior may be an important tool for discriminating between inflationary and defect-based 
models. Furthermore, perturbations generated by scale-dependent variations in the equation 
of state, as considered in this paper, may be a generic property of defect-based cosmological 
scenarios. For global defects, however, the generation of entropy perturbations will be due 
to the distribution of Goldstone, rather than gravitational, radiation. Similar calculations 



have been carried out for textures by Crittenden & Turok [11 



V. CONCLUSION 

In this paper we have analyzed the effects of entropy perturbations in the cosmic string 
scenario. These perturbations arise from spatial variations in the density of cosmic string 
loops and gravitational radiation, leading to a field of approximately Gaussian fluctuations 
in the equation of state of the cosmological fluid. This source of anisotropy has not been 
explicitly considered in previous analyses, because the contribution to the cosmic string 
stress-energy tensor by gravitational waves has not been included. By constructing a crude 
model of this source, we have been able to estimate the effects of the resulting perturbations 
on large scale structure formation and the cosmic microwave background. 

Let us comment now on this analytic model. We claim that the distribution function 
used to describe the variations in the density of gravitational radiation and loops, and 
hence the variations in the equation of state, is a reasonable first approximation for the 
purposes of estimating the relative importance of this effect. There are certain shortcomings 
of this construction, among them the uncertainties in the amplitude of the rms fluctuations 
A, the size of the comoving coherence scale ^, and the dependence of the distribution / on 
wavenumber k, which we have condensed into the parameter a. For large scale perturbations, 
the amplitude and shape of the spectrum does not depend strongly on a or ^, as a scale- 
invariant isocurvature power spectrum results. For what we expect to be a reasonable range 
for the amplitude, A ~ 0.6 — 1, the entropy perturbations make a non-negligible contribution 
to the LSS and large-angle CMB spectrum. In summary, we have crudely modeled the 
essential features; a statistically homogeneous, isotropic distribution of fluctuations in the 
equation of state, varying on length scales above the comoving horizon radius. It remains 
to compare this model with the distribution observed in a realistic cosmic string simulation. 
We plan to carry out these advances in the near future. 
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